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ABSTRACT 

The calculation of the thermal stratification in the superadiabatic layers of stellar models with convective envelopes is a long standing 
problem of stellar astrophysics, and has a major impact on predicted observational properties like radius and effective temperature. 
The Mixing Length Theory, almost universally used to model the superadiabatic convective layers, contains effectively one free pa¬ 
rameter to be calibrated -ttmi- whose value controls the resulting effective temperature. Here we present the first self-consistent stellar 
evolution models calculated by employing the atmospheric temperature stratification, Rosseland opacities, and calibrated variable a ml 
(dependent on effective temperature and surface gravity) from a large suite of three-dimensional radiation hydrodynamics simulations 
of stellar convective envelopes and atmospheres for solar stellar composition (Trampedach et al. 2013). From our calculations (with 
the same composition of the radiation hydrodynamics simulations), we find that the effective temperatures of models with the hydro- 
calibrated variable a mi (that ranges between ~1.6 and ~2.0 in the parameter space covered by the simulations) display only minor 
differences, by at most ~30-50 K, compared to models calculated at constant solar a m i (equal to 1.76, as obtained from the same 
simulations). The depth of the convective regions is essentially the same in both cases. We have also analyzed the role played by the 
hydro-calibrated T{t) relationships in determining the evolution of the model effective temperatures, when compared to alternative 
T{t) relationships often used in stellar model computations. The choice of the T(t) can have a larger impact than the use of a variable 
a ml compared to a constant solar value. We found that the solar semi-empirical T(t) by Vernazza et al. (1981) provides stellar model 
effective temperatures that agree quite well with the results with the hydro-calibrated relationships. 
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1. Introduction 

The almost universally adopted method for calculating superadi¬ 
abatic convective temperature gradients in stellar evolution mod¬ 
els is based on the formali sm p rovided by the so-called Mixing 
Length Theory (MLT - lBohm-VitenselfT958h . This formalism is 
extremely simple; the gas flow is made of columns of upward 
and downward moving convective elements with a characteristic 
size, the same in all dimensions, that cover a fixed mean free path 
before dissolving. All convective elements have the same phys¬ 
ical properties at a given distance from the star centre; upward 
moving elements release their excess heat into the surrounding 
gas, and are replaced at their starting point by the downward 
moving elements, that thermalize with the surrounding matter, 
thus perpetuating the cycle. The MLT is a ‘local’ theory, and 
the evaluation of all relevant physical and chemical quantities 
are based on the local properties of each specific stellar layer, 
regardless of the extension of the whole convective region. 


ent MLT ‘flavours'), so that practically the only free parame¬ 
ter to be calibrated is a m i. Stellar evolution calculations for a 
fixed mass and initial chemical composition but with varying 
a m \, produce evolutionary tracks with different 7j.fr (and radius) 
evolution, whereas evolutionary timescales and luminosities are 
typically unchanged. On the other hand, the prediction of ac¬ 
curate values of T e $ (and radii) by evolutionary models and 
stellar isochrones is paramount, among others, to study colour- 
magnitude diagrams of resolved stellar populations, empirical 
mass-radius relations of eclipsing binary systems, and predict re¬ 
liable integrated spectra (and colours) of unresolved stellar sys¬ 
tems (extragalactic clusters and galaxies). 

In stellar evolution calculations the value of a m i is usually 
calibrated by reproducing the radius of the Sun at the solar age 
with an evolutionary solar model. This solar calibrated a m \ is 
then kept fixed in all evolutionary calculations of stars of dif¬ 
ferent masses and chemical compositions. The exact numeri¬ 
cal value of a m \ varies amongst calculations by different au¬ 
thors because variations of input physics and choices of the 
outer boundary conditions affect the predicted model radii and 
7' c ii values, hence require different or m i values t o mat ch the Sun. 
Regarding the various M LT flavours, iPedersen et al.1 (1 19901) and 
Salar is & Cassisil ((2008) have shown how they provide the same 
T eS evolution, once ar m i is appropriately recalibrated on the Sun. 

As discussed by, e.g., IVan denberg et al.1 (1 19961) and 
Salaris et al. (2002) metal poor red giant branch (RGB) mod¬ 
els calculated with the solar calibrated value of a m \ (hereafter 


Both the mean free path and the characteristic size of the 
convective elements are assumed to be same for all convective 
bubbles, and are assigned the same value A = a m \H p , the so- 
called ‘mixing length’). Here a m \ is a free parameter (assumed 
to be a constant value within the convective regions and along all 
evolutionary phases), and 77 p is the local pressure scale height. 
There are additional free parameters in the MLT, that are gener¬ 
ally fixed a priori (versions of the MLT with different choices 
of these additional parameters will be denoted here as differ- 
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a„,i, o) are able to reproduce the T e « of samples of RGB stars 
in Galactic globular clusters, within the error bars on the es¬ 
timated T eS . On the other hand, several authors find that vari¬ 
ations of a m i with respect to (r ln i 0 are necessary to repro¬ 
duce -just to give some examples- the red edge of the RGB 
in a sample of 38 nearby Galactic disc stars with radii de¬ 
termined from interferometry dPiau et alJl201 ll) . the asteroseis- 
mically constrained radii of a sam ple of main sequence (MS) 
Kepler targets dMathur e t al.||201 2j), observations of binary MS 
stars in the H vades (lYildiz et al.l l2006h and the a Cen system 
(Yildizl 20071). a lo w-mass pre-MS eclipsing binary in Orion 
( Stassun et al.fi2004h . 

The MLT formalism provides only a very simplified de¬ 
scription of convection, and there have been several attempts 
to int r oduce non-localit y in the MLT (see, e.g.. [Grossma n et al.l 
1 19931: IPeng et al.l 120 06, and references therein). These ‘refine¬ 
ments’ of the MLT are often complex and introduce addi- 
tional free parameter s to be calibr ated. The alternative mod el by 
ICanuto & Mazziteliil dl99lh and lCanuto & Mazzitellil d 19921) in¬ 
cludes a spectrum of eddy sizes (rather than the one-sized con¬ 
vective cells of the MLT) and fixes the scale length of the convec¬ 
tive motions to the distance to the closest convective boundary. 
Recently iPasetto et al.l d2014l) have presented a new non-local 
and time-dependent model based on the solution of the Navier- 
Stokes equations for an incompressible perfect fluid, that does 
not contain any free parameteiQ 

An alternative approach to model the superadiabatic layers 
of convective envelopes is based on the computation of real¬ 
istic multidimensional radiation hydrodynamics (RHD) simu¬ 
lations of atmospheres and convective envelopes -where con¬ 
vection emerges from first principles- that cover the range of 
effective temperatures f7Lrr), surface gravities (g), and compo¬ 
sitions typical of stars with surface convection. These simu¬ 
lations have reached nowadays a high level of sophistication 
(see, e.g.. iNordlund et ah] 12009 ) and for ease of implementa¬ 
tion in stellar evolution codes, their results can be used to pro¬ 
vide an ‘effective’hydro-calibration of a m \, even though RHD 
simulations do not confirm the basic MLT picture of columns 
of convective cells. After early attempts from rather crude two- 
dime nsional (2D) and three-dimensional (3D) s imulations (see, 
e.g. JDeupree & Varneill980l : lLvdon et al.lll992l) . a first compre¬ 
hensiv e R HD calibration of a m i was presented by Ludwig et al.l 
am and lFrevtag et al.l d 1 9991) . These authors found from their 
simulations that the calibrated a m i varies as a function of metal- 
licity, g and 7).n . IFrevtag & Sal arTsI (1 1 9991 ) applied this RHD cal¬ 
ibration of a m i (plus T(t) relations computed from the same 
RHD models, and Rosseland opacities consistent with the RHD 
calculations) to metal poor stellar evolution models for Galactic 
globular cluster stars, and found that the resulting isochrones for 
the relevant age range have only small T e g differences along the 
RGB (of the order of ~50 K) with respect to isochrones com¬ 
puted with a solar calibrated value of or m |. 

In the last years a number of grids of 3D hydrodynamics sim¬ 
ulations of surface convection have been published, and from the 
point of view of stellar model calculations it is very im portant to 
study whether the results by iFrevtag & Salarisl d 19991) are con¬ 
firmed or drastically changed. 

iTanner et al.l d2013hl) and iTanner et al.l d2013al) have pre¬ 
sented a grid of simulations employing in the optically thin lay¬ 
ers a 3D Eddington solver (iTanner et al.l 12012T) . Their calcula¬ 
tions cover four metallicities (from Z=0.001 to Z=0.04), but just 


1 This model has not been implemented yet in any stellar evolution 
computations. 


a few T eS values at constant log(g)=4.30for each Z, plus a subset 
of models at varying He mass fraction for Z=0.001 and Z=0.02. 
These authors studied the properties of convection with vary¬ 
ing T eff and chemical composition in these solar-like envelopes. 
ITanner et al.l d2014l) extracted metallicity-dependent T(j) rela¬ 
tions from these same simulations, and employed them to high¬ 
light the critical role these relations play when calibrating a m \ 
with stellar evolution models. 

iMagic et~afl d2013l ) have published a very large grid of 3D 
RHD simulations for a range of chemical compositions. Their 
grid covers a range of T e jf from 4000 to 7000 K in steps of 
500 K, a range of log(g) from 1.5 to 5.0 in steps of 0.5 dex, 
and metallicity, [Fe/H], from -4.0 to +0.5 in steps of 0.5 and 1.0 
dex. These models have been employed bv lMagic et al.1 (}2015h to 
calibrate <+ m | as function of g, T e g and [Fe/H]. They found that 
a m \ depends in a complex way on these three parameters, but 
in general a m \ decreases towards higher effective temperature, 
lower surface gravity and higher metallicity. So far iMagic et al.l 
( 2015 ) have provided only fitting formulae for a m \ but not pub¬ 
licly available prescriptions for the boundary conditions and in¬ 
put physics. _ _ 

Very recently iTrampedach et al.l ( 2013 1) produced a non¬ 
square grid of convective atmosphere/envelope 3D RHD sim¬ 
ulations for the solar chemical composition. The grid spans a 
T e ff range from 4200 to 6900 K for MS stars around log(g)=4.5, 
and from 4300 to 5000 K for red giants with log(g)=2.2, the 
lowest surface gravity available. In ITrampedach et aO (12014bi) 
the horizontal and temporal averages of the 3D simulations were 
then matched to ID hydrostatic equilibrium, spherically sym¬ 
metric e nvelope models to calibra te a m \ as function of g and 
r eff (see Tramneda ch et alJ l2014bl for details about the cali¬ 
bration procedure ). Moreover, the sa me RHD simulations have 
been employed bv lTrampedach et ahl (j2014al) to calculate g- and 
T e ff-dependent 77r) relations from temporal and r (Rosseland 
optic al depth) averaged temp eratures of the atmospheric lay¬ 
ers. fTrampedaeFetlil] (j2014bj) also provide routines to calculate 
their g- and r e ff-dependent RHD-calibrated o- m | together with 
their computed T{ r) relations, and Rosseland opacities consis¬ 
tent with the opacities used in the RHD simulations. This en¬ 
ables stellar evolution calculations where boundary conditions, 
superadiabatic temperature gradient and opacities of the con¬ 
vective envelope are consistent with the RHD simulations. It is 
particularly important to use both the RHD-calibrated a m \ and 
T{t) relations, because the T e g of the stellar evolution calcula- 
tions depends on b oth these inputs (see, i.e.. lSalaris et al.ll2002t 
Tanner et al. 2014, and references therein). 

Thanks to this consistency between RHD simulations and 
publicly available stellar model inputs, we present and dis- 
cuss in this paper the firs t stellar evolution calculations where 
ITrampedach et al.l d2014bh 3D RHD-calibration of is self- 
co nsistently included in the evolutionary code. In the same vein 
as IFrevtag & Salarisl ( 1999 ). we focus on the effect of the cali¬ 
brated variable a-,,,! on the model r e ff, compared to the case of 
calculations with fixed <r m L Q (as determined from the same RHD 
simulations). Self consistency of opacity and boundary condi¬ 
tions is paramount to assess correctly differential effects, given 
that the response of models to variations of a mi depends on their 
Teff, that in turn depends on the absolute value of a m \, opacities 
and boundary conditions. We also address the role played by the 
RHD calibrated T(t ) relations in the determination of the model 
T e g when compared to other widely used relations. 

Section 2 describes briefly the relevant input physics of the 
models, while Sect. 3 presents and compares the resulting evo¬ 
lutionary tracks. A summary and discussion close the paper. 
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2. Input physics 

All stellar evolution calculations presented here have been per- 
forme d with the BaSTI (a BA g of Stellar Tracks and Isochrones) 
code (Pietrinfe rni et al.ll2004l). for a chemical composition with 
Y- 0.245, Z=0.018 and the lGrevesse & Noelsl (119931) metal mix¬ 
ture, consistent with the chemical composition of the RHD simu¬ 
lations. Atomic diffusion was switched off in these calculations, 
and co nvective core overshooti ng was included when appropri¬ 
ate (see lPietrinferni et~aDl2004i for details). 

We employed the Rosseland opacities provided by 
iTrampedach et al.l (1201 4bl) . Th e low temperature opa cities (for 
log(7')<4.5) are very close to iFerguson et al. (2005 b with the 
excep tion of the region with log(7')<3.5, where IFerguson et al.l 
(2005) calculations are higher be cause of the inclusion o f 
the effect of water molecules (see [Tr amnedac het al.l l2014al) . 
Opacity Project (iBadnell et al.l l2005h calculations are used 
for log(r)>4.5. As for the surface boundary conditions^we 
employ ed the T(t) relations computed by Trampedach et al J 
( 2014al) . who provided a routine that calculates, for a given T e g 
and surface gravity, the appropriate generalized Hopf functions 
q(r), related to the T (r) relation by 


qir) = 


4/T(t)\ 2 * 4 
3 \T eff ) 


~ T 


( 1 ) 


In our calculations we fixed to r tl -=2/3 the transition from 
the T(t) integration of the atmospheric layers (with r as inde¬ 
pendent variable) to the integration of the full system of stellar 
struct ure equations. As mentio ned bvlTrampeda ch et alJ ([2014a) 
and ITrampeda ch et al.l (l2014bl) . these RHD-based T(r) relations 
can be employed to model also the convective layers in the opti¬ 
cally thick part of the envelope, together with modified expres¬ 
sions for the temperature gradients in t he superadiabatic regions, 
and an appropriate rescaling of r (see ITrampedach et al.ll2014aL 
for details). In our calculations we have compared the radiative 
(V,ad) and superadiabatic (V -obtained with the RHD calibrated 
values of o^ni) temperature gradients determined along the upper 
part of the convective envelopes from the standard stellar struc¬ 
ture equations and MLT (down to r=100, the upper limit for the 
routine calculating the Hopf functio ns), with V ra d and V calcu - 
lated according to Eqs.35 and 36 of ITrampedach et all ( 2014al) . 
respectively. We have found that in all our models the differ¬ 
ences between these two sets of gradients are much less than 1% 
between r tr and r = 100. 

We have calculated also test models by changing Ttr be¬ 
tween 2/3 and 5 (with the appropriate rescaling of r if convec- 
tion ap pears in the atmosphere integration, see Trampedach et al. 
l2014ai) . and obtained identical evolutionary tracks in each case. 

Regarding the value of a m \, the same routine for the Hopf 
functionsjprovides also the RHD calibrated value of (for the 
iBohm- Vitensel 1958L flavour of the MLT) for a given 73^ and sur¬ 
face gravity, that we employed in our calculations. Uncertainties 
in the cali brated a m i values are of t he order of +0.02 - 0.03 (see 
Table 1 of ITrampedach et ai~j|20 1 4ht) P1. 

The only difference in terms of input physics between the 
atmosphere/envelope RHD calculations and our models is the 
equation of state (EOS). The RHD simulations employed the 
MHD (iDaeppen et alJ fl988l) EOS, that is not the same EOS 

2 We did not include any turbulent pressure in the convective enve¬ 

lope, as the a m i RHD calibration was performed in a way that works for 

standard stellar evolution models without this extra contribution to the 
pressure (R. Trampedach private communication, see also Sect.4 from 
iTrampedach et al.ll2014bl) 


used in the BaSTI calculations (see iPietrinferni et al.ll2004l) . To 
check whether this can cause major differences in the models, 
we have calculated envelope models for the same g-73n pairs of 
the RHD simulations, including the RHD-calibrated a m i, T(t) 
relations, and RHD opacities. We have then compared the re¬ 
sulting depths of the convecti on zones (dry., in units of stellar 
radius) with what obtained by ITrampedach et all ( 2014b ) from 
their RHD-calibrated ID envelope models, that used the same 
input phys ics (including EOS) of th e RHD calculations (see 
Table 1 of ITrampedach et all 12014bl) . We found random (non 
sy stematic) differences of de v by at most just 2-3% compared 
to ITrampedach et all (12014b ) results. 

3. Model comparisons 

As mentioned in the introduction, the RHD simulations cover a 
non-square region in the g-73ff diagram, as displayed in Fig. Q] 
(the region enclosed by thick solid lines), ranging from 4200 to 
6900 K on the MS, and from 4300 to 5000 K for RGB stars with 
log(g)=2.2. The MLT calibration results in an a m \ varying from 
1.6 for the warmest dwarfs, with a thin convective envelope, up 
to 2.05 for the coolest dwarfs in the grid. In between there is a 
triangular plateau of a m \ ~ 1.76, where the Sun is located. The 
RHD simulation for the Sun provides a m \ ~ 1.76+0.03. The top 
panel of Fig. |T| displays the results of our evolutionary model 
calculations in the g-T e g diagram (from the pre-MS to the lower 
RGB), for masses M=0.75, 1.0, 1.4, 2.0 and 3.0M G respectively, 
that cover the full domain of the RHD simulations. 

The thick solid lines denote the reference set of models, cal¬ 
culated with the varying a m i calibration, and the run of o-mi along 
each individual track is displayed in the lower panel. The dotted 
part of each sequence in this panel denotes the a m \ values ex¬ 
trapolated by the calibration routine, when the models are out¬ 
side the region covered by the simulations but still retain a con¬ 
vective region. This happens along the pre-MS evolution of the 
M=O.85M 0 track, and for the 2.0 and 3.0M o calculations along 
the subgiant phase. Apart from the pre-MS stages of the two 
lowest mass models, the evolution of a m i spans a narrow range 
of values, between ~1.6 and ~1.8. 

To study the significance of the variation of a m \ for the 
model T eB , we have calculated evolutionary models for the same 
masses, this time keeping a m i=a m / 0 =1.76 along the whole evo¬ 
lution. The results are also displayed in Fig.Q] 

A comparison of the two sets of tracks clearly shows that the 
effect of a varying a m \ is almost negligible. The largest differ¬ 
ences are of only ~30 K along the RGB phase of the 3.0M o track 
(solar a tracks being hotter because of a higher a m \ value com¬ 
pared to the calibration) and ~50 K at the bottom of the Hayashi 
track of the 1.0M o track (solar a m \ tracks being cooler, because 
of a lower a m i). In all other cases differences are smaller, and 
often equal to almost zero. 

For all stellar masses we found that the mass fraction of He 
dredged to the surface by the first dredge up -that depends on 
the maximum depth of the convective envelope at the beginning 
of the RGB phase- is the same within 0.001, between constant 
a m i' 0 and variable a m \ models. We have then compared the lu¬ 
minosity of the RGB bump -that also depends on the maximum 
depth of the convective envelo pe at the first dredge up (see, e.g., 
ICassisi & Salarisl fl9971 1201 3l and references therein)- for the 
0.75 and 1.0M o models. We found that the luminosity is un¬ 
changed between models with constant and variable a m \. This re¬ 
flect the fact that the depth of the convective envelope is the same 
between the two sets of models throughout the MS phase to the 
RGB, until the end of the calculations. Small differences appear 
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Fig. 1. Stellar evolution tracks in the g-T e g diagram for the la¬ 
belled masses. The region enclosed by the thick black bound¬ 
ary is the g-Teff range covered by the RHD simulations. Thick 
solid lines denote fully consistent calculations with the RHD 
calibrated variable a m \ and 77 r) relationships. The lower panel 
displays the evolution of a m \ along each track. The dotted por¬ 
tion of each sequence denotes the region where the a m \ values 
are extrapolated. Dashed lines in the upper panel display tracks 
calculated with a constant a m \ =a m i,e and the calibrated 77 r) re¬ 
lationships. 


only during the pre-MS. To this purpose, for the 1.0M o models 
we have additionally checked the surface Li abundance that sur¬ 
vives the pre-MS depletion. This provides information about the 
evolution of the lower boundary of the convective envelope dur¬ 
ing this phase, where according to Fig. |TJ a m \ shows the largest 
difference from We found that models with constant a m / i0 
display after the pre-MS a Li abundance just 9% higher than the 
reference results. 

We have then analyzed the role played by the 77r) relation¬ 
ships computed from the RHD calculations in determining the 
model 7’ e ff (the role played by the boundary conditions in deter¬ 
mining the r e ff of stella r models is espec i ally crucial for very - 
low-mass stars, see i.e. lAllard et alJfl997t iBrocato et al.lll998h . 
Figure [2] displays the ratio T/7’ e ff as a function of r as predicted 
by the calibrated relationships for atmospheres/envelopes with 
7^=4500 K and 6000 K (log(g)=3.5). The RHD calibrated 7’(r) 
relationships contain values for the Hopf function that vary with 
r and 7’eff (also with g, to a lesser degree). The variation with 
7^- is obvious from the figure. For these two temperatures the 
largest differences appear at the layers with r between —0.1 and 
~ -1, where stellar model calculations usually fix the transi¬ 
tion from the atmosphere to the interior. The same Figure [2] also 
displays the results for the traditional Eddington approximation 
to the grey atmosphere, and the solar semi-empirical 77r) rela- 
tionships by K r ishna- Swamy (KS - Krish na~SwamvllT966l) and 
IVernazza et alJ (Il98lh -their Model C for the quiet sun, here¬ 
inafter VAL c T(t). In these latter cases the ratio T/7’ e ff does not 
depend on 7’eff. It is easy to notice that around the photospheric 



log(r) 

Fig. 2. Comparison of the ratio T/T e g as a function of the op¬ 
tical depth t as predicted by different T(t) relationships. The 
thick solid and dotted lines display the RHD-calibrated relation¬ 
ships for r e ff=4500 K and 6000 K (log(g)=3.5 in both cases), 
respectively. The thin solid, dash-dotted and dashed lines show 
the ratios obtained from the Eddington, KS and VAL C 7’(r) rela¬ 
tionships, respectively. 


layers the RHD relationships are in between the Eddington grey 
and VAL c T(t). The most discrepant relationship is the KS one. 

Figure [3] displays the results of evolutionary calculations 
with this set of T{t) relationships. We compare here the refer¬ 
ence set of models for 1.0 and 1.4M Q and varying a m \, with con¬ 
stant a m [,Q models calculated using the Eddington, KS and VAL C 
T(t) relationships (we set r„ = 2/3 also for these calculations). 

Differences between these new sets of models at constant 
a m ig) and the reference calculations are larger than the case of 
Fig. HI because of the differences with the RHD 77r) relation¬ 
ships. On the whole, the VAL C T(t) (coupled to the RHD cali¬ 
brated ar m / 0 ) gives the closest match to the self-consistent refer¬ 
ence calculations. The largest differences appear for the 1.0M o 
along the pre-MS; models calculated with the VAL C relation are 
~50 K cooler, approximately the same as the case of a m ; i0 and 
the RHD T(t). For the same l.OM 0 track the T e g differences 
along the RGB are at most equal to 10 K, and at most -40 K 
along the MS. Differences for the 1.4M 0 track are smaller. 

As for the Eddington T(r), the resulting tracks are gener¬ 
ally hotter than the self-consistent RHD-based calculations. The 
largest differences amount to -40-50 K along MS and RGB of 
the two tracks. Comparisons with Fig. Q] show that, from the 
point of view of the resulting model 7’eff, the Eddington T(t) 
differs more -albeit by not much- than the VAL C one from the 
RHD-calibrated relations. 

The worse agreement is found with the KS T(r). For both 
masses the RGB is systematically cooler by -70 K, and the pre- 
MS by -80-100 K, whilst the MS of the l.OM 0 calculations is 
cooler by -120 K, a difference reduced to -50 K along the MS 
of the 1.4M 0 track. 
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Fig. 3. As the upper panel of Fig. |T| but for evolutionary tracks 
with masses equal to 1.0 and 1.4M G respectively. Thick solid 
lines denote fully consistent calculations with the RHD cali¬ 
brated variable a m \ and 7 (r) relationships. Dotted, dash-dotted 
and dashed lines display tracks calculated with constant a m / i0 
and the Eddington, KS and VAL C T(r), respectively. 

In case of the 1.0M o models we checked again the luminos¬ 
ity level of the RGB bump, the surface He mass fraction after the 
first dredge up and the amount of surface Li after the pre-MS de¬ 
pletion. For all these three T (r) relations the evolutionary mod¬ 
els display a RGB bump luminosity within A(L/L G ) <0.01, and 
a post-dredge up He abundance within 0.001 of the fully consis¬ 
tent result. In all these calculations the convective envelope has 
the same depth throughout the MS phase and along the RGB. 
Differences along the pre-MS are highlighted by the amount of 
Li depletion during this phase. We found that models calculated 
with the Eddington T(t) have 9% less Li after the pre-MS, com¬ 
pared to the reference calculations. Comparing this number with 
just the effect of a constant a m i ?Q discussed above, we derive that 
the use of this T{t) decreases the surface Li abundance by -20% 
compared to the use of the Hopf functions determined from the 
RHD simulations. In case of the VAL C T(r), the net effect is to 
increase the surface Li after the pre-MS by -10% compared to 
the RHD T(r), and also the KS relation causes a similar increase 
by -11%. 

3. 1. The standard solar model 

We close our analysis by discussing the implications of the RHD 
results and the choice of the T(t) relation, on the calibration of 
the standard solar model. As well known, in stellar evolution it 
is customary to fix the value of a m ; j0 (and the initial solar He 
and metal mass fractions) by calculating a 1M 0 stellar model 
that matches the solar bolometric luminosity and radius at the 
age of the Sun, with the additional constraint of reproducing 
the present metal to hy drogen mass fraction Z/X ratio (see, e.g., 
iPietrinferni et all2004 for details). The accuracy of the derived 
solar model can then be tested against helioseismic estimates of 


the depth of the convective envelope and the surface He mass 
fraction. It is also well established that solar models without 
microscopic diffusion cannot properly account for some helio¬ 
seismic constraints, hence solar models are routinely calculated 
by including microscop ic diffusion of He and metals (see, e.g., 
iPietrinferni et alj[2004 for a discussion and references). 

We have first calculated a standard solar model (with the 
same input physics and solar metal distribution of the calcula¬ 
tions discussed above) employing both the variable a m \ and the 
T(t) relations from the RHD results. Given that microscopic dif¬ 
fusion decreases with time the surface chemical abundances of 
the model, the initial solar Z (and Y) need to be higher than the 
present one. This means that we had to employ the RHD results 
also for chemical compositions not exactly the same as the com¬ 
position of the RHD simulations. 

We found that it is necessary to rescale the RHD a m \ cali¬ 
bration by a factor of just 1.034 to reproduce the solar radius. 
This implies a m i >0 =1.82, extremely close, within the error, to 
the RHD value ami G = 1.76±0.01 (range )+0.03( calibration uncer- 
tainty) obtained by (Tramne dach et al.ll2014bl) . Given our pre¬ 
vious results, it is also obvious that a solar model calibration 
with fixed ami and the RHD T(t) relations provides the same 
n^mi.0— 1 -82. 

Solar calibrations with the KS, VAL C and Eddington T(r) 
relations have provided a m i i0 =2.11, 1.90 and 1.69, respectively. 

In all these calibrated solar models the initial solar He mass 
fraction (Ti n i j0 - 0.274) and metallicity ((Zj n i i0 - 0.0199) are 
essentially the same, as expected. Also, the model present He 
mass fraction in the envelope (7 0 = 0.244) and the depth of the 
convection zone (dcz = 0.286 R Q ) are the same for all calibra- 
tions and in agreement with the helioseismic values (see, e.g., 
iDziembowski et all 19951 iBasu & Antiall 19971) . 

If we take ami >0 obtained with the RHD T(r) as a reference, 
the lower value obtained with the Eddington T(r), and the larger 
values obtained with both the VAL C and KS relations (in increas¬ 
ing order) are fully consistent with the results of Fig. [3] In that 
figure VAL C and KS T(r) MS models are increasingly hotter than 
the reference RHD calculations (hence increasingly larger a m i 
values are required to match the reference MS) whereas the use 
of the Eddington T(r) produces models cooler than the reference 
MS (hence lower a m i value are needed to match the reference 
MS). 

4. Summary and discussion 

We have presented the first self-consistent stellar evolu- 
tion calculations that em ploy the variable a m i and 74 r) by 
iTrampedach et alj (l2014bl) . based on their 3D RHD simulations. 
Our set of evolutionary tracks for different masses and the same 
chemical composition (plus consistent Rosseland opacities) of 
the RHD simulations, cover approximately the entire g-7 01 pa¬ 
rameter space of the 3D atmosphere/envelope calculations. 

We found that, from the point of view of the predicted 7’ e ff 
(plus the depth of the convective envelopes and amount of pre- 
MS Li depletion), models calculated with constant RHD cali¬ 
brated a m i. Q = 1.76 are very close to, and often indistinguish¬ 
able from, the models with variable a m \. Maximum differences 
are at most -30-50 K. Th is result is similar to the conclusions by 
iFrevtag & Salarisl dl999t) . based on 2D RHD simulations at low 
metallicities. 

At first sight this may appear surprising, given that the full 
range of a m i spanned by the RHD calibration is between -1.6 
and -2.0. However, one has to take into account the following 
points: 
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1. The derivative A7' e ff/Aa' m i for stellar evolution models de¬ 
pends on the absolute value of a m \, and decreases when ar m i 

increases). 

2. The variation of AT e g for a given Aa m i depends also on the 

mass extension of the convective and superadiabatic regions. 

It is therefore clear that the decrease of with increasing 
r eff cannot have a major effect because of the thin convective 
(and superadiabatic) layers of models crossing this region of the 
g-Teff diagram. The large variations An,,,! ~0.2 (see lower panel 
of Fig. Q]) along the lower Hayashi track of the 1.0M o model 
is also not very significant (~50 K) because of the decreased 
extension of the surface convection and the reduced AT e g/Aa^ 
at higher a m \. 

It is however very important to remark here that the detailed 
structure of the superadiabatic convective regions is not suitably 
reproduced either by a m / i0 or by a variable ir m \, and that the full 
results from RHD models need to be employed whenever a de¬ 
tailed description of the properties of these layers is needed. 

To some degree the role played by the RHD T(t ) is more 
significant. Constant 0 = 1-76 stellar evolution models be¬ 
come systematically cooler by up to ~ 100 K along the MS, 
pre-MS and RGB when the widely used KS T(j) relation is 
used, compared to the self consistent RHD-calibrated calcula¬ 
tions. Evolutionary tracks obtained employing the Eddington 
and VAL C T(t) relationships are much less discrepant, and stay 
within ~50 K of the self-consistent calculations. Even from the 
point of view of pre-MS Li depletion the VAL C T( r) causes only 
minor differences, of the order of 10%. A similar difference is 
found with the KS relation, whilst a larger effect of ~2()% is 
found with the Edd ington T(r). _ 

An extension of lTrampedach et aD ( 2013 ) simulations to dif¬ 
ferent metallicities is necessary to extend this study and test the 
significance of the a m \ variability (and Hopf functions)^over a 
larger parameter space. The 3D RH D simulations by Magic et al J 
(2013) and the a m \ calibration by iMagic et al J (l2015h cover a 
large metallicity range, and predict an increase of <r ln | with de¬ 
creasing metal content, but at the moment it is not possible to 
properly include this calibration in stellar evolution calculations, 
for the lack of available T(t) relations and input physics consis¬ 
tent with the RHD simulations. At solar metallicity the general 
beha viour of ami with g and TVff seems to be qualitatively simi¬ 
lar to lTramnedach et aD (12014hl) calibration. However, the range 
of a m i i s shifted to higher values compared to Tramnedach et alJ 
(1201 4bl ) results, due probably to different input physics and the 
different adopted solar chemical composition. 

Acknowledgements. We are grateful to R. Trampedach for clarifications about 
his results, and for making available the routine to calculate Rosseland opacities. 
We also thank J. Christensen-Dalsgaard for interesting discussions on this topic, 
and the anonymous referee for comments that have improved the presentation of 
our results. SC warmly acknowledges financial support from PRIN-INAF2014 
(PI: S. Cassisi) 


Daeppen, W., Mihalas, D., Hummer, D. G., & Mihalas, B. W. 1988, ApJ, 332, 
261 

Deng, L„ Xiong, D. R.. & Chan, K. L. 2006, ApJ, 643, 426 
Deupree, R. G. & Varner, T. M. 1980, ApJ, 237, 558 

Dziembowski, W. A., Goode, P. R.. Pamyatnykh, A. A., & Sienkiewicz, R. 1995, 
ApJ, 445, 509 

Ferguson, J. W., Alexander, D. R., Allard, F.. et al. 2005, ApJ, 623, 585 
Freytag, B., Ludwig. H.-G., & Steffen, M. 1999, in Astronomical Society 
of the Pacific Conference Series, Vol. 173, Stellar Structure: Theory and 
Test of Connective Energy Transport, ed. A. Gimenez, E. F. Guinan, & 
B. Montesinos, 225 

Freytag, B. & Salaris, M. 1999, ApJ, 513, L49 
Grevesse, N. & Noels, A. 1993, Physica Scripta Volume T, 47, 133 
Grossman, S. A., Narayan, R., & Arnett, D. 1993, ApJ, 407, 284 
Krishna Swamy, K. S. 1966, ApJ, 145. 174 
Ludwig, H.-G., Freytag, B„ & Steffen, M. 1999, A&A, 346, 111 
Lydon, T. J., Fox, P. A., & Sofia, S. 1992, ApJ, 397, 701 
Magic, Z„ Collet, R„ Asplund, M„ et al. 2013, A&A, 557, A26 
Magic, Z„ Weiss, A., & Asplund, M. 2015, A&A, 573, A89 
Mathur, S., Metcalfe. T. S., Woitaszek, M., et al. 2012, ApJ, 749, 152 
Nordlund, A., Stein, R. F., & Asplund, M. 2009, Living Reviews in Solar 
Physics, 6, 2 

Pasetto, S., Chiosi, C„ Cropper, M„ & Grebel, E. K. 2014, MNRAS, 445, 3592 
Pedersen, B. B., Vandenberg, D. A., & Irwin, A. W. 1990, ApJ, 352, 279 
Piau, L„ Kervella, P„ Dib, S., & Hauschildt, P. 2011, A&A, 526, A100 
Pietrinferni, A., Cassisi, S., Salaris, M., & Castelli, F. 2004, ApJ, 612, 168 
Salaris, M. & Cassisi, S. 2008, A&A, 487, 1075 
Salaris, M„ Cassisi, S., & Weiss, A. 2002, PASP, 114, 375 
Stassun, K. G., Mathieu, R. D., Vaz, L. P. R.. Stroud, N., & Vrba, F. J. 2004, 
ApJS, 151, 357 

Tanner, J. D., Basu, S., & Demarque, P. 2012, ApJ, 759, 120 
Tanner, J. D., Basu, S., & Demarque, P. 2013a, ApJ, 778, 117 
Tanner, J. D., Basu, S., & Demarque, P. 2013b, ApJ, 767, 78 
Tanner, J. D., Basu, S., & Demarque, P. 2014, ApJ, 785, L13 
Trampedach, R., Asplund, M., Collet, R., Nordlund, A., & Stein, R. F. 2013, 
ApJ, 769, 18 

Trampedach, R., Stein, R. F., Christensen-Dalsgaard, J., Nordlund, A., & 
Asplund, M. 2014a, MNRAS, 442, 805 

Trampedach, R., Stein, R. F., Christensen-Dalsgaard, J., Nordlund, A., & 
Asplund, M. 2014b, MNRAS, 445, 4366 
Vandenberg, D. A., Bolte, M., & Stetson, P. B. 1996, ARA&A, 34, 461 
Vemazza, J. E., Avrett, E. H., & Loeser, R. 1981, ApJS, 45, 635 
Ytldtz, M. 2007, MNRAS, 374, 1264 

Ytldtz, M„ Yakut, K„ Baki§, H„ & Noels, A. 2006, MNRAS, 368, 1941 


References 

Allard, F., Hauschildt, P. H., Alexander, D. R.. & Starrfield, S. 1997, ARA&A, 
35, 137 

Badnell, N. R„ Bautista, M. A., Butler, K„ et al. 2005, MNRAS, 360, 458 
Basu, S. & Anda, H. M. 1997, MNRAS, 287, 189 
Bohm-Vitense, E. 1958, ZAp, 46, 108 

Brocato, E., Cassisi, S., & Castellani, V. 1998, MNRAS, 295, 711 
Canuto, V. M. & Mazzitelli. I. 1991, ApJ, 370, 295 
Canuto, V. M. & Mazzitelli. I. 1992, ApJ, 389, 724 
Cassisi, S. & Salaris, M. 1997, MNRAS, 285, 593 

Cassisi, S. & Salaris, M. 2013, Old Stellar Populations: How to Study the Fossil 
Record of Galaxy Formation 


6 





















